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Abstract 

NuMI neutrino beam is constructed to aim at the MINOS detector in 
Soudan mine. Neutrinos emitted at angles 10 — 20 mrad with respect to 
the beam axis create an intense beam with a well defined energy, dependent 
on the angle. Additional surface detectors positioned at the transverse dis- 
tance of several kilometers from the mine offer an opportunity for very precise 
mesurements of the neutrino oscillation parameters. The mixing matrix ele- 
ment \U&\ can be measured down to a value of 0.0025 with the exposure of 
the order of 20 kton x years. 
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Introduction 



A discovery of the neutrino oscillations by the SuperK|I| experiment has changed 
the neutrino physics. The existence of the phenomenon is no longer in question and 
the focus of the experiments has changed from the 'discovery' to the 'elucidation'. 
The next round of experiments, running or under construction, will provide an in- 
dependent verification of the SuperK result and it is likely confirm that neutrinos 
oscillate mostly into v T with neutrinos with a characteristic frequency corresponding 
to Am 2 ~ 1.5 -4 x 10~ 3 eV 2 . 

Questions related to neutrino masses and oscillations are of a fundamental na- 
ture therefore a significant effort is dedicated to a possible design of the future 
facilities: high intensity conventional neutrino beams, often referred to as 'super- 
beams', and/or a novel facility, 'neutrino factory', where a very intense mixed beam 
of z/^'s and u e 's is produced from decay of muons circulating in a storage ring. 

In this note we examine a potential extension of the utilization of the NuMI 
beam, being constructed at Fermilab, by addition of new detectors. These new 
detectors could offer a major enhancement of the physics reach of the MINOS de- 
tector, especially in conjunction with the anticipated improvements of the Fermilab's 
accelerator complex. 

2 The NuMI Beams 

Neutrino beam constructed at Fermilab has a flexible design 0. Focusing elements, 
consisting of two parabolic horns, can be be moved with respect to the target and 
to each other thus focusing different momentum bytes of the produced pions and 
giving rise to the neutrino beams of different energies. 

The resulting beams are presented in Fig.Q showing the event spectra for a nom- 
inal 1 kton x year exposure at the distance of 730 km from the target. In this paper 
we assume the proton flux on target to be 3.8 x 10 20 per year. 

The horn acceptance is defined in terms of the relative momentum spread 
producing the neutrino beam with the energy spread increasing with the beam 
energy. The long tail of the high energy neutrinos in the low and medium energy 
beams is due to the bare target component of the beam created by high momentum 
pions traversing the opening in the magnetic horns. 
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Figure 1: High, medium and low energy options for NuMI beams. Histograms 
represent the rates of CC events expected at the distance of 730 km in a detector 
with 1 kton fiducial volume. 

3 The Principle of Off-axis Neutrino Beams 

Most of the neutrinos in a conventional neutrino beam are produced in two-body 
decays 7r + — > fi + + v^. A smaller contribution ti the beam is due to decays, 
K + — > /i + + u^, and to the decays of the muons, /i + — > e + X^> e . 

Energy and the flux of the neutrinos produced in two-body decays is uniquely 
determined by the decay angle 9d ec - For small angles they can be expressed as: 
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1 + 7 2 ^, 
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where 

• 7 = — is the Lorentz boost factor of a pion 

• 9 dec is decay angle, i.e. the angle between the pion and the produced neutrino 
directions 

• A is the area of the detector and z is the distance between the decay point 
and the detector. 




Figs.0 and [3] illustrate the relationship between the neutrino fluxes, energies and 
the decay angles for the pions in the energy ranges relevant to the NuMI beam. 
Neutrino flux is maximal in the forward direction, 9 dec — 0, and the neutrino energy 
there is proportional to the energy of the parent pion, E v = 0A3E n . This is the 
reason for placing the neutrino detector 'on-axis'. 

Fig.|3| shows that the sharp peak of the neutrino flux in the forward direction 
changes into a broad maximum for pions at low energies, E n < 5 GeV. Low energy 
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Figure 3: Variation the neutrino flux with the decay angle 6 dec for neutrinos produced 
by 7T 's of different energies. 

neutrino beam, E v < 2 — 3 GeV, is very broad and the detector location is not very 
critical for the energy and intensity of the neutrino flux. This fact makes low energy 
neutrino oscillation experiments relatively insensitive to the details of understanding 
of the underlying neutrino beam. 

The same fact enables a construction of improved, with higher intensity and 
smaller energy spread, low energy neutrino beams. This technique, pioneered in 
the Brookhaven neutrino oscillation experiment proposal[[|], consists of placing a 
neutrino detector at some angle with respect to the conventional neutrino beam. 

Such a detector records approximately the same flux of low energy neutrinos, 
as the one positioned 'on-axis', originating from the decays of low energy pions. 
In addition, though, an 'off-axis' detector records an additional contribution of low 
energy neutrinos from the decay of higher energy pions decaying at the finite angle 

dec- 
All additional attractive feature of the neutrino flux observed at the 'off-axis' 
detector is a kinematical suppression of high energy neutrino component. Detectors 
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placed at different angles with respect to the neutrino beam direction are therefore 
exposed to an intense narrow-band neutrino with the energy defined by the detector 
position. 

In contrast with the traditional narrow-band neutrino beams the off-axis detec- 
tors enable a simultaneous experiments with beams at different energies/locations. 
At the same time the 'off-axis' neutrino flux is much higher than the one achievable 
with the traditional narrow-band beam. 

4 Off-axis Beams at NuMI 

Neutrino beam described in the previous section is directed towards a MINOS de- 
tector located in the Soudan mine. This beam, especially in its low- and medium- 
energy version is very broad and its energy spectrum changes significantly with the 
position/angle of the neutrino detector. 

We consider a hypothetical situation, where the 'on-axis' MINOS detector is 
complemented by three additional detectors positioned at the transverse distance 
of 5, 10 and 20 km away from the central detector. For a sake of an example we 
examine a case of 10 kton x year exposure. 

4.1 Low Energy Beam Case 

Fig.fl shows the expected (in 'no oscillations') case rates of z/ M CC events in the four 
detectors. 

As expected, the neutrino spectrum is much narrower at the 'off-axis' detectors 
than the one at the central detector, with the peak position shifting to the lower 
energies with the increasing distance from the beam axis. The high(er) energy tail 
of the spectrum is due to high momentum pions produced at relatively large angles 
and decaying shortly after the production target and due to kaon decays, as shown 
in Fig.|. 

4.2 Medium Energy Beam Case 

Expected neutrino events spectra at the 'off-axis' detectors in the case of the medium 
energy NuMI beam are shown in Fig.^[ 

This configuration seems to be close to the optimal running conditions. The 
neutrino flux in the 3 — 4 GeV region at the distance of 5 km is higher than the 
peak of the low energy beam 'on-axis', with the high energy tail greatly reduced. At 
the distance of 10 km the flux of the neutrinos in 1 — 2 GeV region is much higher, 
whereas the contribution of higher energy neutrinos is suppressed in comparison 
with the same detector exposed to the low energy beam in 'on-axis' position. 
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Figure 4: Low energy beam option: neutrino event spectra for 10 kton*years expo- 
sure at the detectors on-axis and at 5,10 and 20 km distance. 



4.3 High Energy Beam Case 

High energy beam is derived from high energy pions produced in the target. At 
these high energies a sharp drop of the neutrino flux intensity, related primarily to 
the large factor 7 = E^/m^ in the EqPJ leads to the effective beam intensities at 
the 'off-axis' detectors significantly reduced in comparison with the low and medium 
energy beams, as shown in FigTj] 

5 Disappearance Experiment with the Off-axis 
Detectors 

A disappearance experiment consisting of several detectors positioned at different 
distances from the neutrino beam axis has several attractive features: 

• average energy of the neutrino beam can be selected by the position of the 
detectors. A suitable choice of the detectors location will enable detection 
of a oscillation pattern whereby the observed event rate goes down and up 



S 




Figure 5: Contribution of pion- and kaon-produced neutrinos to the neutrino flux 
at the detectors on-axis and at 5,10 and 20 km distance, low energy beam. 

with the increasing distance of the detector from the beam and hence with the 
decreasing neutrino energy, as illustrated on Fig.[|. 

• measurement of the disappearance of CC events can be easily extended to 
very low energies, E v < 1 GeV . At these low energies the CC events will be 
primarily of quasi-elastic type, + n — > pT + p and will consist of a single 
penetrating muon track. In case of a wide band beam in the 'on-axis' detector 
such events will be contaminated by single pion events produced by the NC 
interactions of neutrinos at higher energies. 

• a judicious choice of the detector location in conjunction with the small energy 
spread of the main peak of the neutrino flux enables a very precise determi- 
nation of the mixing angle sm 2 2#23- As an example, in case of Am| 3 = 
3 x 10~ 3 GeV a 10 kton detector positioned at the distance of 10 km and 
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Figure 6: Medium energy beam option: neutrino event spectra for 10 kton*years 
exposure at the detectors on-axis and at 5,10 and 20 km distance. 



exposed to the neutrino beam produced by 8 x 10 protons on target will 
observe 94 CC events instead of the expected 821 events. The number of 
'missing' events provides a direct determination of the mixing angle sin 2 2023 
with the statistical accuracy Asm 2 20 2 3 < 0.0014. 

• the average energy of the neutrinos in the peak of the spectrum is determined 
from two-body kinematics and the detector position, independent of the de- 
tector resolution. This is a significant advantage over the wide-band neutrino 
beam, where the kinematical effect of the masses of the produced particles and 
nuclear effects introduce systematic difference between the observed (even with 
an ideal detector) total energy of the final state and the energy of the incoming 
neutrino. This effect is likely to limit the ultimate precision of the determina- 
tion of Am| 3 by the neutrino oscillation experiments, baring a breakthrough 
in the understanding of the intra-nuclear cascades and other nuclear effects. 
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6 Systematics I: Understanding of the Beam Spec- 
tra 

Well defined beam energy, determined completely by the detector position and a 
small energy spread make the 'off-axis' beam potentially a very powerful tool for 
studies of neutrino oscillations. For this purpose it is necessary, however, to be able 
to predict the expected neutrino beam, its spectrum and the overall normalization 
with adequate precision. 

The far detector spectra in the 'on-axis' can be determined from the event spectra 
observed at the near detector. Does the experiment with 'off-axis' detector require 
a dedicated near experiments near the neutrino source? Such an experiment would 
be complicated by the extended nature of the neutrino source. 

Fortunately, such an experiment is not necessary. The neutrinos observed in 
the near detector and in any of the far detectors are derived from the same parent 
beam, formed by the focusing elements. This common origin implies strong corre- 
lations between the beams observed at different locations. The correlation matrix 
is primarily determined by the geometry of the neutrino beam and by the details 
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Figure 8: Expected CC event rates it the on- and off- axis detectors for 10 
kton*years exposure in the medium energy beam. The histogram represents 'no- 
oscillation' case. The cross-hatched distribution shows the expected rate in case of 
oscillations with full mixing and Am 2 = 3.0 x 10~ 3 eV 2 . 



of the focusing system, thus permitting an unambiguous prediction of the neutrino 
flux at the arbitrary position using the spectrum observed at the near detector [|]. 
A residual dependence on the hadron production model is very small, both for the 
case of the low and medium energy beams, as indicated by the Figs.0 and [TO. 



7 Systematics II: Backgrounds 

The principal limitation of a sensitivity of a possible v e appearance experiment 
will be related to the v e and V e component of the beam. They are produced in 
K + — ► 7T°e + u e and fi + — ► e + v e 17^ decays. The angular distribution of these three- 
body decays is somewhat different from the decays producing the main component 
of the beam. 

A relative fluxes of the i/„ and v e neutrinos at different detector positions are 



show in Fig.|Tl| for the low energy beam case and in Fig.[LJ for the medium energy 
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Figure 9: Dependence of the predicted neutrino spectrum on the hadron production 
model for the case of the low energy beam. Four overlayed histograms represent the 
expected spectra at different detector positions derived from the spectrum observed 
at the near detector using different hadron production models. 



beam. 

The v e flux in the 'off-axis' detectors is of the order of 0.5% of the z/ M flux and 
this background is somewhat reduced in comparison with the 'on-axis' detector. 



8 v e Appearance Experiment 

The primary goal of the v e appearance experiments is a measurement of the 
element of the MNS neutrino mixing matrix or the mixing angle sinO^. The ex- 
periment is expected to detect a small number of u e CC interactions among the 
neutrino interactions registered at the far detector. In the following we assume that 
the experiment is carried out at such a combination of L and E v that the oscillations 
generated by Am| 3 dominate and that the value of the Am^ is measured precisely 
in the disappearance experiment. A number of the v e interactions is proportional 
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Figure 10: Dependence of the predicted neutrino spectrum on the hadron produc- 
tion model for the case of the medium energy beam. Four overlayed histograms 
represent the expected spectra at different detector positions derived from the spec- 
trum observed at the near detector using different hadron production models. 



to the neutrino mixing parameters via 



/rIN 1 97Am 2 T 

dE^sm 2 23 = sin 2 29 etl N^ ear (3) 

where 



sin 2 26 eix = -sin 2 26 13 = 2 \U e3 \ 2 (4) 

The v e CC interactions are characterized by an absence of a long, penetrating 
muon track and by a presence of an electromagnetic energy cluster due to the out- 
going electron. There are several sources of backgrounds yielding events with very 
similar characteristics: 
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Figure 11: and v e induced CC event rates at the on- and off-axis detectors, low 
energy beam. The upper histogram represents the rates, the lower one shows the 
v e rates. 



• high-y CC interactions with a significant fraction of the hadronic energy 
produced in a form of 7r°'s. 

• NC interactions with a significant fraction of the hadronic energy produced in 
a form of 7r°'s 

• CC interaction of z/ r 's (resulting from — > v T oscillations) with a subsequent 
decay r — > e. This background is particularly troublesome if the mixing angle 
sin 2 28 efl is much smaller than sin 2 26 T ^. 

An additional, irreducible, contribution to the sample of potential v e CC inter- 
actions is produced by an intrinsic v e component of the neutrino beam. 

A sensitivity of the measurement of the mixing parameters will be limited by 
the statistical fluctuations of the background event sample, even if the overall nor- 
malization of the background is known. It is, therefore, very important to minimize 
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Figure 12: and v e induced CC event rates at the on- and off-axis detectors, 
medium energy beam. The upper histogram represents the rates, the lower one 
shows the v e rates. 



the size of the background sample, while the sample of detected u e CC interactions 
due to the oscillations is maximized. 

The optimization procedure involves several factors: 

• a judicious choice of the beam energy and spectrum. An optimal neutrino 
beam would have a large flux concentrated around E opt = ^— — Such a 
beam would have the following effects: 

— a number of v^s which have 'oscillated' away, jV"* sappear , in Eq|| would 
be maximized 

— background due to high-y CC interactions would be minimized 

— background from the NC interactions would be minimized if the neu- 
trino spectrum does not extend beyond the region where the oscillation 
probability is very large 
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a choice of the distance of the detector from the neutrino source. To eliminate 
the background due to CC u T interactions it is desirable that the optimal 
energy E opt is below the kinematical threshold for r production. 

a suitable choice of the analysis cuts. Oscillations-induced v e events will be 
concentrated in a narrow energy range around E opt . A large fraction of NC 
events will have an observed energy much lower than the neutrino energy due 
to the escaping neutrino, hence a requirement that the detected energy of the 
event exceeds a certain fraction of the E opt may significantly reduce the NC 
background. The intrinsic u e component of the beam tends to have much 



harder energy spectrum (see Figs. 11 and 12), hence a requirement that the 



observable energy does not exceed a certain value will significantly reduce this 
background. 




2 

Figure 13: Sensitivity of the appearance experiments to \U e3 \ min for low energy beam 
exposure 20 kt*years. 

A potential of the v e appearance experiment with the 'off-axis' detectors is shown 
in Figs. [13] and [14] for the case of low and medium energy beams respectively. 
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2 

Figure 14: Sensitivity of the appearance experiments to \U e3 \ min for medium energy 
beam exposure of 20 kt*years. 

2 

Sensitivity of the experiment to the | U e s \ is taken to be equal to the 3 a fluctuation of 
the background. It is assumed that the ir° rejection power of the neutrino detector is 
sufficient to reduce the NC background to the level below the intrinsic u e component 
of the beam. The size of the exposure is taken to correspond to 20 kton*years at 
the 100% electron identification efficiency 

A detector situated at the distance of 10 km from the beam axis offers the best 
opportunity for the measurement of the |?7 e 3| 2 . Although a potential reach of the 
detector located at the distance of 5 km is somewhat better for larger values of Am| 2 , 
the NC background would be much larger there, thus necessitating a substantially 
better detector to ensure the adequate background rejection. 

An example of the actual event rates for the case of medium energy beam and 
oscillations corresponding to Am^ = 0.003 eV 2 is shown in Table |l|. 
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Table 1: Examples of v e appearance experiments for Am^ = O.OOSeV 2 . E min — E max 
denote the optimal energy range, no osc./osc. are the expected numbers of CC 
interactions without and with oscillations, \Uez\ min corresponds to the value of the 
mixing parameters producing an effect equal to 3 a fluctuation of the background. 



Det. position, R 


E m in E max , GeV 


no osc. 


Ve bckg. 


v n osc. 


\U \ 2 

1 e 3 1 rain 


km 


1.0 - 7.5 


13460 


60.0 


9380 


0.0030 


5 km 


1.0 - 4.5 


5420 


24.0 


2480 


0.0025 


10 km 


1.5 - 2.5 


1380 


4.4 


80 


0.0024 


15 km 


1.0 - 1.5 


320 


1.0 


100 


0.0068 


20 km 


0.5 - 1.5 


200 


1.2 


120 


0.0206 



9 Neutrino Detectors Issues 

The NuMI neutrino beam is directed towards the MINOS detector located in the 
Soudan mine. Construction of underground facilities at the off-axis locations is not 
practical, therefore the design of the detectors and their locations must be opti- 
mized to enable precise measurements despite the presence of cosmic rays induced 
backgrounds. 

Transverse distance of the detector from the beam axis determines the average 
energy of the neutrino beam. Thus the optimal detector position depends on the 
value of the Am 2 and on the nature of the desired measurement. 

For example, if Am 2 = 3 x 10~ 3 eV 2 then a detector at the transverse distance of 

10 km from Soudan would be located at the maximum of the oscillation and it would 
be an optimal one for the measurement of the mixing angle and for the detection 
of v e appearance. An additional detector at the distance of 20 km could be used to 
demonstrate the oscillation pattern. 

Optimization of a possible detector location may include the variation of the 
distance of the putative detectors from Fermilab. Sensitivity of the neutrino oscil- 
lations appearance experiments does not depend on the experiment baseline as long 
as the beam energy is much higher than the one corresponding to the oscillation 
maximum, as reduction of the flux is compensated by the increased oscillation 

probability (-§) 2 - 

The situation is different when the combination of the experiment baseline L 
and the beam energy E corresponds to the oscillation maximum. The change of the 
sensitivity of the appearance experiment with the reduction of the baseline of the 
experiment is a combination of an increase of the neutrino flux as and a reduction 
of the oscillation probability, if the detector is located at the same angle with the 
respect to the neutrino axis or a change of the neutrino flux with the increased angle 
if the detector is located at a larger angle to lower the neutrino energy to correspond 
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to the oscillation maximum at the new location. This interplay between the detector 
longitudinal and transverse positions is illustrated in Fig.[L5| 
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Figure 15: Rates of CC event rates in the detectors located at different locations. 
Rates correspond to the medium energy beam, 10 kton x year exposure 

The NuMI neutrino beam has an interesting property of passing, at the distance 
of about 600 km, under the large body of water, the Lake Superior. It allows for 
an unique neutrino oscillation experiment utilizing the same detector exposed to 
the narrow-band neutrino beam of different energies. Such an experiment mounted 
on a barge floating the the Lake Superior can be moved to an arbitrary position, 
thus permitting the optimization of the neutrino beam as the understanding of 
the neutrino oscillation progresses. An experiment can be also built as submerged 
module(s) towed by the barges to the requested positions. Such a scenario can 
significantly reduce the cosmic muon background. 
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